Non-zero macroscopic magnetization in half-metallic antiferromagnets 

at finite temperatures 
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Combining density functional theory calculations with many-body Green function technique, we 
reveal that the macroscopic magnetization in half-metallic antiferromagnets does not vanish at finite 
temperature as for the T = limit. This anomalous behavior stems from the inequivalent magnetic 
sublattices which lead to different intra-sublattice exchange interactions. As a consequence, the 
spin fluctuations suppress the magnetic order of the sublattices in a different way leading to a 
ferrimagnetic state at finite temperatures. Computational results are presented for the half-metallic 
antiferromagnetic CrMnZ (Z = P, As, Sb) semi-Heusler compounds. 
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Half-metallic antiferromagnets (HM-AFMs) are con- 
sidered to be the most promising class of materials for 
spintronics applications.^'^ A HM-AFM material is not 
antiferromagnetic in the usual sense of the term; it is 
a special case of a ferrimagnet with compensated sub- 
lattice magnetization. The existence of the gap in one 
of the spin channels (either up or down) leads to the 
complete cancellation of the magnetic moments at zero 
temperature with a 100% spin polarization of the charge 
carriers at the Fermi level. In conventional AFMs the 
macroscopic spin polarization is zero due to the spin ro- 
tational symmetry with the exception of the compounds 
with broken inversion symmetry like a— Fe2 03 in which 
spin-orbit gives rise to weak ferromagnetism (0.002 /i^) 
due to the canting of the magnetic momentsi^ The HM- 
AFM materials provide several advantages in device ap- 
plications with respect to the half-metallic ferromagnets 
(HM-FMs). For example, they would be perfectly sta- 
ble spin-polarized electrodes in a junction device. These 
materials do not give rise to stray fields, and thus no 
magnetic domain walls are formed. Besides applications 
in spintronics, HM-AFMs provide a possibility of "single 
spin superconductivity" due to the spin triplet {S = 1) 
pairing in metallic channeli^ 

The possible existence of HM-AFM was pointed out 
by van Leuken and de Groot in 1995i^ Based on first- 
principles calculations authors proposed CrMnSb and 
VrMnFcsSbyln as candidates for HM-AFMs. Later, 
Pickett suggested that also the cubic double perovskites 
La2VCu06, La2MnV06, and La2MnCo06 are HM- 
AFMs.^ Since then substantial effort has been de- 
voted to find materials with HM antiferromagnetic 
characteristics. Predicted promising systems include 
double perovskites^ii^, thiospinels^, tetrahedrally co- 
ordinated transition-metal-based chalcopyriteaii, fuU- 
Heusler alloys'^^!-'^ and monolayer superlatticesJ^ Not 
only ordered but also several disordered systems have 
been shown to be HM-AFMs J^ii^^J^ii^i^o In spite of 
substantial theoretical efforts for designing materials 
with HM-AFM characteristics and the study of their 
ground state electronic and magnetic properties, only 



few works exist addressing the exchange interactions 
and magnetic phase transition temperatures in these 
systems] ^^1^° 

In this Communication we investigate the temperature 
dependence of the magnetization in HM-AFMs employ- 
ing the many-body Green function technique^ ^ within 
Tyablikov decoupling schemei^ For computational pur- 
poses we consider CrMnZ (Z = P, As, Sb) semi-Heusler 
compounds which are the simplest systems among the 
predicted HM-AFMs with two magnetic atoms per unit 
cell and which are compatible with the existing semicon- 
ductors technology. However, present findings are valid 
for more complicated systems like double perovskites^ 
or diluted antiferromagnetic semiconductors i^^^^S We 
show that the macroscopic magnetization in these ma- 
terials does not vanish at finite temperature in contrast 
to the zero temperature limit and conventional AFMs. 
This peculiar behavior originates from the inequivalent 
magnetic sublattices in HM-AFMs which lead to differ- 
ent intra-sublattice exchange interactions and, as a conse- 
quence, spin fluctuations suppress the magnetic order of 
the sublattices in a different way. Thus, at finite temper- 
ature sublattice magnetizations do not compensate each 
other and all three compounds show ferrimagnetic behav- 
ior which seems to be contradictory to our knowledge on 
finite temperature properties of the magnetic materials. 
However, this seemingly counterintuitive results can be 
explained by an analysis of the electronic structure and 
exchange interactions in these systems. 

Ground state calculations are carried out using the 
augmented spherical wave method within the general- 
ized gradient approximation to the exchange correlation 
potential. Details of the computational scheme can be 
found in Ref. [H. To provide the basis for further consid- 
erations we start with a brief discussion of the electronic 
structure of the CrMnZ compounds. Like several Heusler 
alloys these systems have theoretical equilibrium lattice 
constants (see Table IJ) close to the ones of the zincblende 
semiconductors (GaP, GaAs). Thus, calculational results 
will be presented for the latter case since these semi- 
conductors might be considered as possible substrates 
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TABLE I: Lattice parameters and spin magnetic moments (in 
fis) for CrMnZ (Z — P, As, Sb). The calculated equilibrium 
lattice constants are 5.44 A, 5.71 A and 6.08 A for Z = P, As 
and Sb, respectively, close to the experimental lattice param- 
eters of GaP, GaAs and InP. 



Compound 


a(A) 


mcr 


niMn 


mz 


mcell 


CrMnP 


5.45[GaP] 


L80 


-L83 


0.03 


0.00 


CrMnAs 


5. 65 [GaAs] 


2.52 


-2.54 


0.02 


0.00 


CrMnSb 


5.87[lnP] 


2.74 


-2.76 


0.02 


0.00 



to grow these materials All compounds under study 
have 18 valance electrons per unit cell and calculations 
show that the total spin moments given in Table U are 
exactly zero in agreement with the Slater- Pauling be- 
havior for ideal half-metals.^^ Simultaneously, the spin 
magnetic moments of Cr and Mn atoms are antiparallel 
and these compounds are ferrimagnets with compensated 
sublattice magnetizations. The origin of the HM gap in 
Heusler alloys has been well understood and the reader 
is referred to Ref. [2^ since the same discussion is valid 
for the present systems. 

To study interatomic exchange interactions we map 
the complex itinerant electron problem onto a classical 
Heisenberg Hamiltonian H = ~ Er,r' -^rr'SrSr' 
where /iR ^ i^R' and the indices /i and ly represent dif- 
ferent sublattices. R and R' are the lattice vectors speci- 
fying the atoms within the sublattices, and is the unit 
vector pointing in the direction of the magnetic moment 
at site (/i,R). Heisenberg exchange parameters J^^, are 
calculated employing the frozen-magnon technique as de- 
scribed in Ref. [23. Extensive investigations on the multi- 
sublattice Heusler alloys have shown that there are sev- 
eral exchange interactions which coexist and which are 
mixed together j^i^iSli^ To simplify the discussion let us 
write the total magnetic exchange field acting on the sub- 
lattice n as Jl'^^^^ - 4(;2t + •Al^'d'^Lct + jlndLct where the 
first two terms represent the direct and indirect exchange 
coupling between different sublattices. The former (di- 
rect coupling) provides the leading contribution to the 
total exchange coupling and determines the character of 
the magnetic state while the latter contributes little and 
is not so important. On the other hand, the last term, 
i.e., the intra-sublattice coupling is of particular impor- 
tance because it is responsible for the appearance of a 
net macroscopic magnetization at finite temperature in 
HM-AFMs. It should be noted here that in reahty the 
situation is not so simple, the exchange field acting on 
the sublattices should be determined from the solution 
of a matrix equation. 

Let us start with the discussion of the last term. 
Because of the large distance between Cr-Cr (Mn-Mn) 
atoms this coupling is indirect mediated by the conduc- 
tion electrons. In Fig. [1] we present the calculated Cr- 
Cr and Mn-Mn Heisenberg exchange parameters as a 
function of the distance. As seen, due to inequivalent 
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FIG. 1: (Color onhne) Intra-sublattice Cr-Cr and Mn-Mn 
exchange interactions in CrMnZ (Z — P, As, Sb) as a function 
of distance. In the insets we show spin-resolved total density 
of states around Fermi level. 



magnetic sublattices the pattern of Cr-Cr and Mn-Mn 
exchange parameters show very different behavior. Al- 
though the former has ferromagnetic character, the latter 
shows oscillatory behavior so that the ferromagnetic and 
antiferromagnetic contributions partly compensate each 
other giving rise to a small contribution into the total 
exchange field acting on the Mn sublattice. In Table |TT] 
we present the intra-sublattice ( = J2'r=£o Jqr) ^^^^ 
as the inter-sublattice {Jq'^ = X^j^ •^qr) on-site exchange 
coupling parameters. The on-site Cr-Cr and Mn-Mn ex- 
change couplings are rather different and this difference 
will be refiected as a net macroscopic magnetization at 
finite temperatures (see Fig. [J). It is worth to note that 
in conventional AFMs this coupling is the same for both 
sublattices. The increase of the strength of the exchange 
interactions and correspondingly of the critical temper- 
atures in the P-As-Sb sequence can be explained by the 
increase of the magnetic moments (see Table |l|. More- 
over, due to the presence of the HM gap the exchange 
interactions quickly decayi^i^l 

In contrast to the intra-sublattice exchange interac- 
tions, the inter-sublattice (Cr-Mn) ones behave very dif- 
ferently. Due to the smaller Cr-Mn distance a very 
strong antiferromagnetic nearest-neighbor direct inter- 
action takes place between these atoms which is about 
one order of magnitude larger than the nearest-neighbor 
Cr-Cr coupling and is responsible for the formation of 
the ferrimagnetic state. The interactions between fur- 
ther nearest neighbors are very small. The on-site Cr- 
Mn exchange couplings presented in Table |TT] are about 
4 and 20 times larger than the Cr-Cr and Mn-Mn 
ones, respectively. Note that the ferromagnetic intra- 
sublattice exchange interactions further stabilize the fer- 
rimagnetic order. The antiferromagnetic coupling of the 
Cr and Mn atoms can be qualitatively explained on 
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TABLE II; On-site exchange parameters (in mRy) and esti- 
mated critical temperatures (in K) within RPA for quantum 
(t[?') and classical (t[?') spins for CrMnZ (Z = P, As, Sb). 



Compound 


■,[Cr-Crl 


-.[Mn-Mn] 

•'o 


■[[Cr-Mn] 

•'o 


T»'(K) 


t[?'(K) 


CrMnP 


3.80 


0.03 


-12.99 


2530 


1264 


CrMnAs 


5.55 


0.84 


-19.50 


2610 


1566 


CrMnSb 


5.96 


1.06 


-20.16 


2986 


1792 



the basis of the following facts: First, half- filled shells 
tend to yield a strong trend towards antiferromagnetism 
and second, exchange coupling in 3d transition met- 
als obeys the semi-phenomenological Bethe-Slater-Neel 
curve which predicts antiferromagnetism in the case of 
small interatomic distances.^- Indeed, the Cr-Mn dis- 
tance ((i[cr-Mn] = 2.36 — 2.54A) in the present systems is 
comparable with the Cr-Cr distance (d[cr-Cr] = 2. 52 A) in 
the antiferromagnetic bcc Cr and both magnetic atoms 
posses half-filled 3c? shells, thus, antiferromagnetic Cr- 
Mn coupling is expected. 

With calculated exchange parameters in hand, now 
we can study temperature dependence of the magneti- 
zation employing the methods of statistical mechanics 
to the Heisenberg Hamiltonian. We use the many-body 
Green function technique^^ within Tyablikov decoupling 
scheme^ [also known as the random-phase approxima- 
tion (RPA)] as described in Ref. [U. Note that RPA 
takes into account only transverse spin fluctuations and 
the spin-flip Stoner excitations (longitudinal spin fluctu- 
ations) are neglected. However, available experimental 
data on Heusler alloys have shown that these latter exci- 
tations are well separated in energy from the former one 
(spin waves) due to large exchange splitting A (A ~ 2 — 3 
eV).^^ In addition to this the presence of HM gap pre- 
vents spin-flip transitions. Thus, Stoner excitations do 
not not play an important role in thermodynamics of the 
present systems and the RPA method is well grounded. 
We consider both classical-spin and quantum-spin cases. 
In the classical-spin calculations the obtained values of 
the magnetic moments (see Table H]) are used while for 
the quantum mechanical case we assign integer values to 
the atomic moments: 2/zb {S = 1) (for Cr and Mn) in 
CrMnP and 3ms {S = 3/2) in CrMnAs and CrMnSb. In 
quantum-spin case the thermal average of the sublattice 
magnetization is given by the Callen's expression'^*' 

^ (S, <&,)(! + <1>J^^"+^ + (^;. + 1 + <&,)< 

(l + $^)2S,+r_($^)2S, + l 

where is an auxiliary function. As the <I>p depends 
on the (sp as well as (sj) to be determined, the above 
equation forms a self-consistency problem to be solved 
by iteration. Note that for classical spins the Callen's ex- 
pression is reduced to the Langevin function.^ In Fig. [2] 
we present the calculated temperature dependence of the 
sublattice and total magnetization for CrMnZ (Z = P, 
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FIG. 2: (Color online) The calculated temperature depen- 
dence of the sublattice and total magnetization for CrMnZ 
(Z = P, As, Sb). The calculations are performed for both clas- 
sical and quantum Hamiltonians. The temperature is given 
in reduced form and the total magnetization is scaled up by 
a factor of 5. 

As, Sb). The temperature is given in reduced form and 
the total magnetization is scaled up by a factor of 5. As 
seen, the spin fluctuations suppress the magnetic order 
of the sublattices differently, i.e., the magnetization of 
the Mn sublattice decreases faster than the Cr one and 
as a result the total magnetization does not vanish at 
flnite temperature in contrast to the T = limit. This 
behavior can be traced back to the different exchange 
flelds acting on the sublattices (see Table |lll . The to- 
tal magnetization shows non- monotonous behavior i.e., 
flrst, it increases with increasing temperature up to 0.8 
Tc and then decreases much faster and flnally becomes 
zero at the critical point. For quantum mechanical case 
the calculated total magnetic moment around 0.8 Tc 
is 0.1 - 0.16/is (0.006 - O.Ol/is around room tempera- 
ture) which is about two orders of magnitude larger than 
the spin-orbit coupling induced weak magnetic moment 
(0.002 [Ib) in a— Fe203.— We should also note that in 
these systems not only the ideal case of zero macroscopic 
magnetization but also half-metallicity is limited to the 

T = 0,2 -32.33 

The nature of the spin (quantum or classical) plays an 
important role in the temperature behavior of the mag- 
netization curves. In the quantum case, the magnetiza- 
tion drops slower than the classical case and thus, the 
calculated Tc values are larger by a factor of [S + 1)/S' 
entering the RPA expression (see Ref. [sil ) . Note that this 
factor becomes unity for classical spins ((5 -I- l)/5 ^ 1 
for S — > cxd). Another important point which is outside 
the scope of the present work is that in both treatments 
we use the exchange parameters estimated within the 
picture of classical atomic moments. However, it is pos- 
sible that the values of the exchange parameters must 
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be modified for the use in the quantum-mechanical cal- 
culations. In general, the classical calculation provides 
reasonable values of the critical temperature compared 
with experiment while the quantum mechanical treat- 
ment gives better form of the temperature dependence 
of the magnetization.— The calculated critical tempera- 
tures within RPA are presented in Table [TTl We notice 
that the predicted Tc values of the CrMnAs and CrMnSb 
are even higher than the fee Co which possesses the high- 
est critical temperature (1400 K) among all known mag- 
netic materials. This is not surprising, because available 
experimental and theoretical data have shown that the 
critical temperatures of HM ferromagnets (or ferrimag- 
nets) scales linearly with the average value of the mag- 
netic moment per atom in the unit cell.^^ In this respect 
Co2FeSi possesses the largest average magnetic moment 
per atom of 1.5/iB with an experimental Tc of 1100 K. 
However, in CrMnAs (CrMnSb) the average value of the 
absolute magnetic moment per atom is 1.68/^^ (l.SS/is) 
which is larger than the corresponding value in Co2FeSi, 
and thus such high critical temperatures are expected. 
Finally we should note that so far discussion is based 
on the assumption that CrMnZ compounds possess Clb- 
type ordered crystal structure and the effect of disor- 
der is completely ignored. However, in reality disorder 
exist in various forms like defects, anti-sites or atomic 



swaps which reduce not only the spin polarization at 
the Fermi level but also the magnetic phase transition 
temperature.— For example migration of the Cr atoms 
to the vacant sublattice is expected to reduce Tc sub- 
stantially since as shown in Rcf. 13 the L2i-type Cr2MnZ 
compounds have Tc values around room temperature. 

In summary, we have studied the temperature depen- 
dence of the magnetization in HM-AFM CrMnZ (Z = P, 
As, Sb) compounds employing the many-body Green 
function technique. We have shown that the macroscopic 
magnetization in these systems does not vanish at fi- 
nite temperature in contrast to the T = limit. This 
anomalous behavior stems from the incquivalcnt mag- 
netic sublattices in HM-AFMs which lead to different 
intra-sublattice exchange interactions and, as a result, 
spin fluctuations suppress the magnetic order of the sub- 
lattices in a different way. Thus, at finite temperatures, 
the sublattice magnetizations do not compensate each 
other and all three compounds show ferrimagnetic behav- 
ior. Moreover, the combination of large HM gaps, high 
Tc values and very small macroscopic magnetization 
around room temperature makes CrMnAs and CrMnSb 
promising candidates for spintronics applications. 

Fruitful discussions with L. Sandratskii, I. Galanakis 
and Ph. Mavropoulos are greatly acknowledged. 



* Electronic address: "e.sasioglu®fz-juelich.de" 

^ Warren E. Pickett and Jagadeesh S. Moodera, Physics To- 
day 54, 39 (2001). 

^ M. I. Katsnelson, V. Yu. Irkhin, L. Chioncel, and A. I. 
Lichtenstein, and R.A. de Groot, Rev. Mod. Phys. 80, 315 
(2008). 

^ L. M. Sandratskii and J. Kiibler, Europhys. Lett. 33, 447 
(1996). 

* V. V. Mazurenko, V. I. Anisimov, Phys. Rev. B 71, 184434 

(2005) . 

^ Warren E. Pickett, Phys. Rev. Lett. 77, 3185 (1996). 
® H. van Leuken and R. A. de Groot, Phys. Rev. Lett. 74, 
1171 (1995). 

Warren E. Pickett, Phys. Rev. B 57, 10613 (1998). 

* J. H. Park, S. K. Kwon, and B. I. Min, Phys. Rev. B 65, 
174401 (2002). 

® Y. K. Wang and G. Y. Guo, Phys. Rev. B 73, 064424 

(2006) . 

1° M. S. Park, S. K. Kwon, and B. I. Min, Phys. Rev. B 64, 

100403(R) (2001). 
" M. Nakao, Phys. Rev. B 74, 172404 (2006). 

S. Wurmehl, H. C. Kandpal, G. H. Fecher, and C. Felser, 

J. Phys.: Condens. Matter 18, 6171 (2006). 

I. Galanakis, K. Ozdogan, E. §a§ioglu, and B. Akta§, Phys. 

Rev. B 75, 172405 (2007). 

M. Nakao, Phys. Rev. B 77, 134414 (2008). 

I. Galanakis, K. Ozdogan, E. §a§iogIu, and B. Akta§, Phys. 

Rev. B 75, 092407 (2007). 

D. Kodderitzsch, W. Hergert, Z. Szotek, and W. M. Tem- 
merman, Phys. Rev. B 68, 125114 (2003). 
^"^ Yung-mau Nie and Xiao Hu, Phys. Rev. Lett. 100, 117203 



(2008). 

^® H. Akai and M. Ogura, Phys. Rev. Lett. 97, 026401 (2006). 
L. Bergqvist and P. H. Dederichs, J. Phys.: Condens. Mat- 
ter 19, 216220 (2007). 

M. Ogura, C. Takahashi, and H. Akai, J. Phys.: Condens. 
Matter 19, 365226 (2007). 

P. Frobrich and P. J. Kuntz, Physics Reports, 432, 223 
(2006) 

S. V. Tyablikov, Methods of Quantum Theory of Mag- 
netism (Plenum Press, New York, 1967). 
E. §a§ioglu, L. M. Sandratskii, and P. Bruno, Phys. Rev. 
B 70, 024427 (2004). 

Among the considered compounds only CrMnSb is experi- 
mentally synthesized with a crystal structure different than 
the Clf, type semi-Heusler one [J. H. Wijngaard et ai, 
Phys. Rev. B 45, 5395 (1992)]. However, state-of-the-art 
experimental methods for the synthesis of materials such 
as the molecular-beam epitaxy make possible the growth 
of these compounds in metastable phases as thin films or 
multilayers where the substrate determines the lattice pa- 
rameter. 

■^^ I. Galanakis and Ph. Mavropoulos, J. Phys.: Condens. 

Matter 19, 315213 (2007). 
^® J. Rusz, L. Bergqvist, J. Kudrnovsky, and I. Turek, Phys. 

Rev. B 73, 214412 (2006). 

J. Kudrnovsky, V. Drchal, L Turek, and P. Weinberger, 
Phys. Rev. B 78, 054441 (2008). 

Ralph Skomski, Simple Models of Magnetism, (Oxford 
University Press, 2008). 
■^^ P. J. Webster and K. R. A. Ziebeck, in Alloys and Com- 
pounds of d-Elements with Main Group Elements, Part 2, 



5 



edited by H. R. J. Wijn, Landolt-Bornstein, New Series, 

Group III, Vol. 19, part. C (Springer, Berlin, 1988). 

H. B. Callen, Phys. Rev. 130, 890 (1963). 

E. §a§ioglu, L. M. Sandratskii, P. Bruno and I. Galanakis, 

Phys. Rev. B 72, 184415 (2005). 

P. A. Dowben and R. Skomski, J. Appl. Phys. 93, 7948 
(2003). 



M. Lezaic, Ph. Mavropoulos, J. Enkovaara, G. Bihlmayer, 
and S. BIgel Phys. Rev. Lett. 97, 026404 (2006). 
J. Kiibler, G. H. Fecher and C. Felser, Phys. Rev. B 76, 
024414 (2007). 

B. Ailing, S. Shallcross, and I. A. Abrikosov, Phys. Rev. B 
73, 064418 (2006). 



